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1 Introduction 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) is 
one of the two remaining operating hadron colliders, the other being the LHC at CERN. 
Unlike the LHC, which is buried in a deep underground tunnel, RHIC is built in an enclosure 
on the surface which is covered by an earth berm for shielding which can be seen from outer 
space (Fig. [I]). 



Figure 1: NASA infra-red photo of Long Island and the New York Metro Region from space. RHIC is the 
white circle to the left of the word BNL. Manhattan Island in New York City, ^100 km west of BNL, is also 
clearly visible on the left side of the photo. 

*Research supported by U. S. Department of Energy, DE-AC02-98CH10886. 
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BNL is a multipurpose laboratory, quite different in scope from Fermilab and CERN, 
with many “cutting edge” major research facilities in addition to RHIC. Figure [2] shows the 
two newest facilities: the National Synchrotron Light Source II (NSLS II) to come on-line 
October 1, 2014; and the Long Island Solar Farm which has an experimental section as 
well as supplying 32MW peak power to nearby homes in partnership with the local electric 
company. 



Figure 2: Aerial view of BNL with NSLS II, RHIC and the Solar Farm indicated. 

2 News from BNL since ISSP2013 

Although Fiscal Year 2014 started on October 1, 2013 with the U. S. Government shut down 
for the first 16 days due to the lack of an approved budget, the rest of the FY2014 turned out 
very well for BNL and RHIC. At the administrative level, Professor Robert Tribble of Texas 
A&M University, well known for both his physics research and leadership of two important 
U. S. Department of Energy (DoE) Panels in 2005 and 2013, was appointed Deputy Director 
for Science & Technology, effective February 2014. Not long after that, in March, the DoE 
issued a Request For Proposals for a “management and operating (M&O) contractor” for 
BNL, which is owned by the U.S. Government, but run by an M&O contractor. The present 
contractor, BSA, is a partnership of Battelle Memorial Institute, a private non-profit science 
and technology development company, headquartered in Columbus, Ohio, and Stony Brook 
University. BSA has been the M&O contractor at BNL for the past 15 years (out of the BNL’s 
67 year existence), with the “engagement” of six of the world’s leading research universities 
(Columbia, Cornell, Harvard, MIT, Princeton and Yale) who were among the universities 
that formed the founding M&O contractor, Associated Universities Incorporated, along with 
Johns Hopkins, and the Universities of Pennsylvania and Rochester. The new contract starts 
on January 1, 2015, preceded by a maximum 2 month transition phase-in period so should 
be awarded near or soon after November 1, 2014 [lj. 
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The U. S. High Energy Physics bureaucracy was not idle during this period, with the 
release of the “Particle Physics Project Prioritization Panel (P5)” Report to the High energy 
Physics Advisory Panel (HEPAP) on May 21, 2014. The charge of the panel was “to develop 
an updated strategic plan for U.S. high energy physics that can be executed over a 10 year 
timescale, in the context of a 20 year global vision for the held.” Their reasonable top priority 
for constrained budget scenarios was to “Use the Higgs boson as a new tool for discovery” 
which is good news for the U. S. HEP groups working at the LHC at CERN; but lots of 
internal U. S. activities were “redirected” [2j. For the unconstrained budget scenario, all 
they could come up with was: 

• Develop a greatly expanded accelerator R&D program that would emphasize the ability 
to build very high-energy accelerators beyond the High-Luminosity LHC (HL-LHC) 
and ILC at dramatically lower cost. 

• Play a world-leading role in the ILC experimental program and provide critical exper¬ 
tise and components to the accelerator, should this exciting scientific opportunity be 
realized in Japan. 

which IMHO lacks the imagination and drive of previous generations of U.S High Energy 
Physicists who had proposed and were constructing a 40 TeV p+p collider for completion in 
1995 if not for ... [3] . 

Not to be outdone, the new Long Range Planning exercise for U. S. Nuclear Physics was 
initiated in April 2014. 

3 RHIC Operations and accelerator future plans 

Since beginning operation in the year 2000, RHIC, which can collide any species with 
any other species including polarized protons, has provided collisions at 14 different val¬ 
ues of nucleon-nucleon c.m. energy, and ten different species combinations including 

Au+Au, d+Au, Cu+Cu, Cu+Au, U+U, and in 2014 He 3 +Au, if differently polarized pro¬ 
tons are counted as different species. The performance history of RHIC with A+A and 


RHIC Run 

I Year 

Species 

Energy 

Ldt 

Run-1 

2000 

Au+Au 

130 GeV 

1 nb-1 

Run-2 

2001-2 

Au+Au 

200 GeV 

24 ab-1 

Run-2 


Au+Au 

19 GeV 

0.4 pb-l 



p+p 

200 Gev 

150 nb-1 

Run-3 

2002/3 

d+Au 

200 GeV 

2.74 nb-1 



P+P(L) 

200 GeV 

0.35 nb-1 

Run-4 

2003/4 

Au+Au 

200 GeV 

241 ab-1 



Au+Au 

62.4 GeV 

9 ab-1 

Run-5 

2005 

Cu+Cu 

200 GeV 

3 nb-1 



Cu+Cu 

62.4 GeV 

0.19 nb-1 



Cu+Cu 

22.4 GeV 

2.7 ab-1 



P+P(L) 

200 GeV 

3.8 pb-l 

Run-6 

2006 

p+p (T+L) 

200 GeV 

10.7 pb-l 



p+p 

62.4 GeV 

100 nb-1 

Run-7 

2007 

Au+Au 

200 GeV 

813 ab-1 

Run-8 

2007/2008 

d+Au 

200 GeV 

80 nb-1 



P+P (T) 

200 GeV 

5.2 pb-l 



Au+Au 

9.2 GeV 


Run-9 

2009 

P+P (L) 

200 GeV 

16 pb-l 



P+P (L) 

500 GeV 

14 pb-l 

Run-10 

2010 

Au+Au 

200 GeV 

1.3 nb-1 



Au+Au 

62.4 GeV 

100 ab-1 



Au+Au 

39 GeV 

40 ab-1 



Au+Au 

7.7 GeV 

260 mb-1 

Run-11 

2011 

P+P (L) 

500 GeV 

27 pb-l 



Au+Au 

200 GeV 

915 ab-1 



Au+Au 

27 GeV 

5.2 ab-1 



Au+Au 

19.6 GeV 

13.7 M events 

Run-12 

2012 

P+P (L) 

200 GeV 

9.2 pb-l 



P+P (L) 

510 GeV 

30 pb-l 



U+U 

193 GeV 

171 ab-1 



Cu+Au 

200 GeV 

4.96 nb-1 

Run-13 

2013 

P+P (L) 

510 GeV 

156 pb-l 

Run-14 

2014 

Au+Au 

15 GeV 

44.2 pb-l 



Au+Au 

200 GeV 

2.56 nb-1 



He3+Au 

200 GeV 

134 nb-1 


Years 

Beam Species and 
Energies 

Science Goals 

New Systems 
Commissioned i 

2014 

15 GeV Au+Au 

200 GeV Au+Au 

Heavy flavor flow, energy loss, 
thermalization, etc. 

Quarkonium studies 

Electron lenses 

56 MHz SRF 

STAR HFT, MTD 

2015-16 

p+p at 200 GeV 
p+Au, d+Au at 200 

GeV 

Extract r|/s(T) + constrain initial 
quantum fluctuations 

More heavy flavor studies 

PHENIX MPC-EX 
Coherent e-cooling test 

2017 

No Run 


Low energy e-cooling 
upgrade 

2018-19 

5-20 GeV Au+Au 
(BES-II) 

Search for QCD critical point 
and onset of deconfinement 

STAR iTPC upgrade 
Partial commissioning 
of sPHENIX (in 2019) 

2020 

No Run 


Complete sPHENIX 
installation 

STAR forward upgrades 

2021-22 

Long 200 GeV Au+Au 
with upgraded 
detectors p+p, p/d+Au 
at 200 GeV 

Jet, di-jet, y-jet probes of parton 
transport and energy loss 
mechanism, color screening for 
different quarkonia 

sPHENIX 

2023-24 

No Runs 


Transition to eRHIC 


Figure 3: a)(left) Year, species and proton polarization (Longitudinal or Transverse), ■ s / r sJJ, integrated 
luminosity of RHIC runs, b) (right) Future run schedule and new equipment. 
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polarized p+p collisions is shown in Fig. [5^t; and in Fig. [5]b, the plans for future runs are 
shown. 

For this year’s run (2014) the full 3 dimensional cooling including electron lenses for 
partial compensation of the beam-beam tune shift and 56 MhZ storage r.f. for stronger 
longitudinal focusing were implemented which led to a higher initial luminosity and much 
longer lifetime of the beam with a more level luminosity load due to the 3d stochatic cooling 
(Fig. §. The luminosity performance of RHIC with A+A and polarized p+p collisions is 
shown in Fig. [5| Notably, the Au+Au f C dt in 2014 exceeds all previous Au+Au runs 
combined as did the p+p f C dt in 2013. 



Figure 4: Run-14 luminosity vs. storage time compared to Run-7, courtesy Wolfram Fischer. 




0 2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20 

Time [weeks in physics] Time [weeks in physics] 


Figure 5: a) (left) Au+Au performance, where the nucleon-pair luminosity is defined as Ann = A x B x L, 
where L is the luminosity and A , B are the number of nucleons in the colliding species, b) (right) Polarized 
p+p performance. Courtesy Wolfram Fischer. 
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The major future plan for accelerators in Nuclear Physics concerns an electron-ion col¬ 
lider, which if located at BNL will be called eRHIC. A new highly innovative and cost- 
effective design of eRHIC was proposed this year based on a Fixed Focus Alternating Gra¬ 
dient (FFAG) electron accelerator and an Energy Recovery Linac (Fig. [6]). 



• 80% polarized e 

► E = 6.6-21.2 GeV 

• 70% polarized protons 

► E = 25 -275 GeV 

• Ions (d ... U) 

► 10 -110 GeV/u 

• Vs = 30- 145 GeV 

• L® 1-3x10 33 ernes' 1 


• Full use of existing RHIC complex including tunnel & cryo 

• 1.32 GeV Energy Recovery Linac with 99.5% recovery efficiency 

• Novel FFAG lattice allows 16 beam re-circulations using only two beam transport rings 

• Permanent magnet technology is used for the FFAG beamline 

• Initial cost estimate: FY14$: 750M (not including detectors) 


Figure 6: New BNL design for eRHIC with annotations. 


4 Detector issues in A+A compared to p+p collisions 

A main concern of experimental design in RHI collisions is the huge multiplicity in A+A 
central collisions compared to p+p collisions. A schematic drawing of a collision of two 
relativistic Au nuclei is shown in Fig. [7^. In the center of mass system of the nucleus- 



Figure 7: a) (left) Schematic of collision in the N-N c.m. system of two Lorentz contracted nuclei with 
radius R and impact parameter b. The curve with the ordinate labeled da/dn c h represents the relative 
probability of charged particle multiplicity n c h which is directly proportional to the number of participating 
nucleons, 7V part . b) (right) E T distribution in Au+Au at +s^n = 200 GeV from PHENIX [4]. 
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nucleus collision, the two Lorentz-contracted nuclei of radius R approach each other with 
impact parameter b. In the region of overlap, the “participating” nucleons interact with 
each other, while in the non-overlap region, the “spectator” nucleons simply continue on 
their original trajectories and can be measured in Zero Degree Calorimeters (ZDC), in fixed 
target experiments, so that the number of spectators can be measured from which the number 
of participants (A part ) can be determined for symmetric A+A collisions. The degree of 
overlap is called the centrality of the collision, with b ~ 0 , being the most central and 
b ~ 2 R, the most peripheral. The maximum time of overlap is r Q = 2 Rj+c where 7 is the 
Lorentz factor and c is the speed of light in vacuum. The energy of the inelastic collision is 
predominantly dissipated by multiple particle production, where N^, the number of charged 
particles produced, or E T , the energy emitted transverse to the beam direction, is directly 
proportional jS] to N paxt as sketched on Fig. [7k. Thus, N ch and Ex in central Au+Au 
collisions are roughly A times larger than in a p+p collision, as shown in actual events from 
the STAR and PHENIX detectors at RHIC (Fig. [ 8 ]). 



Figure 8: a) (left) A p+p collision in the STAR detector viewed along the collision axis; b) (center) Au+Au 
central collision at yTvT = 200 GeV in STAR; c) (right) Au+Au central collision at s/+R = 200 GeV in 
PHENIX. 


At colliders, the impact parameter b can not be measured directly because charged spec¬ 
tators are swept away from zero degrees by the collider magnets. Instead, the centrality of 
a collision is defined in terms of the upper percentile e.g. top 10%-ile, upper 10—20%-ile, 
of IVch or E T distributions as in Fig. 0>- Unfortunately the “upper” and “-ile” are usually 
not mentioned which sometimes confuses the uninitiated. Also a model is required to derive 
Apart from the measurement so that the derived value of A par t at a collider or the number 
of binary nucleon-nucleon collisions (A co u) is model dependent and may have biases. 

Since it is a huge task to reconstruct the momenta and identity of all the particles pro¬ 
duced in these events, the initial detectors at RHIC | 6 j concentrated on the measurement of 
single-particle or multi-particle inclusive variables to analyze RHI collisions, with inspiration 
from the CERN ISR which emphasized those techniques before the era of jet reconstruction 
(see, for example, Refs. |7] and ®)- There are at present two major detectors in operation 
at RHIC, STAR and PHENIX, and there were also two smaller detectors, BRAHMS and 
PHOBOS, which have completed their program. As may be surmised from Figs. [ 8 ^,b and 
UK STAR, which emphasizes hadron physics, is most like a conventional general purpose 
collider detector, a Time Projection Chamber to detect all charged particles over the full 
azimuth (A0 = 27r) and ±1 units of pseudo-rapidity ( 77 ); while PHENIX (Figs. [ 8 ^ and [ 9 ) 3 ), 
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Figure 9: Actual STAR (a) and PHENIX (b) detectors, compare with Figs[8]D,c. The direction of the beam 
is along the axis of the STAR solenoid; and in (b) between the two spectrometer arms (perpendicular to the 
caption). 

is a very high granularity high resolution special purpose detector: a two-arm spectrometer 
at mid-rapidity, with each arm covering solid angle \rj\ < 0.35, A</> = 90°, together with two 
full-azimuth muon detectors at forward and backward rapidity (1.1 < |// < 2.3) Q For the 
present runs, both STAR and PHENIX have excellent particle identification (PID) capability 
with electromagnetic calorimeters (EMcal) for photon and electron detection and Time of 
Flight for charged hadrons. PHENIX has a Ring Imaging CHerenkov counter for enhanced 
electron detection and triggering and small but full azimuth EM calorimeters (MPC) just 
before each muon arm covering 3.1 < \r]\ < 3.7, while STAR obtains enhanced hadron iden¬ 
tification using dE/dx in the TPC. For the 2014 run, both PHENIX (VTX, FVTX) and 
STAR (HFT) are equipped with micro-vertex detectors for tagging Heavy-Flavor c and b 
quarks via displaced vertices. 

The main objectives of buliding RHIC were i) to discover the Quark Gluon Plasma 
(QGP), which was achieved as I have discussed in detail in review articles based on previ¬ 
ous ISSP proceedings Eli; ii) to measure its properties, which were much different than 
expected, namely a “perfect fluid” of quarks and gluons with their color charges exposed 
rather than a gas. The latest measurements from RHIC continue to be very interesting. 

5 7V c h, E t distributions and constituent-quarks as the 
fundamental elements of particle production 

The first experiment specifically designed to measure the dependence of the charged particle 
multiplicity in high energy p+A collisions as a function of the nuclear size was performed 
by Wit Busza and collaborators at Fermilab using beams of ~ 50 — 200 GeV/c hadrons 
colliding with various fixed nuclear targets. They found the extraordinary result ra that 
the average charged particle multiplicity (N c h)hA in hadron+nucleus (h+A) interactions was 
not simply proportional to the number of collisions (absorption-mean-free-paths), (X co ii) = A 
but increased much more slowly, proportional to the number of participants (N pait ) . Thus, 
relative to h+p collisions (Fig. [lOh) [llj: 

Ra = (Ach )hA / (A^ch )hp = (Apart )hA / (-Apart)fy p = (1 + v)/2 . (1) 

1 The detector is so non-conventional that it made the cover of Physics Today, October 2003. 
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Figure 10: a) R A = {N c h) hA / {N c h ) h as a function of the average thickness of each nucleus given in terms 
of the mean free path, V = (N co n) [11 for 50 and 100 GeV/c h+A collisions; b) Charged particle multiplicity 
density, dN c h/drj, as a function of A (represented by V) for 200 GeV/c p+A collisions [T2l . 


Since the different projectiles, h = 7r + , K + ,p in Fig. 10 1 have different mean free paths, the 
fit to the same straight line in terms of V is convincing. 

The other striking observation (Fig. [ToJd) [T2j was that a relativistic incident proton could 
pass through e.g. v = 4 absorption-mean-free-paths of a target nucleus and emerge from 
the other side; and furthermore there was no intra-nuclear cascade of produced particles (a 
stark difference from what would happen to the same proton in a macroscopic 4 mean-free- 
path hadron calorimeter). In the forward fragmentation region of 200 GeV/c p+A collisions, 
within 1 unit of rapidity from the beam, y heam = 6.06, there was essentially no change in 
dN c h/dr) as a function of A, while at mid-rapidity (?/ ~ y™ = 3.03), dN^/drj increased 
with A together with a small backward shift of the peak of the distribution resulting in 
a huge relative increase of multiplicity in the target fragmentation region, rj < 1 in the 
laboratory system. These striking features of the ~ 200 GeV/c fixed target hadron-nucleus 
data (y/Sjvjv ~ 19.4 GeV) showed the importance of taking into account the time and distance 
scales of the soft multi-particle production process including quantum mechanical effects. 


5.1 The Wounded Nucleon Model 


The observations in Fig. [10] had clearly shown that the target nucleus was rather transparent 
so that a relativistic incident nucleon could make many successive collisions while passing 
through the nucleus, and emerge intact. Immediately after a relativistic nucleon interacts 
inside a nucleus, the only thing that can happen consistent with relativity and quantum 
mechanics is for it to become an excited nucleon with roughly the same energy and reduced 
longitudinal momentum and rapidity. It remains in that state inside the nucleus because 
the uncertainty principle and time dilation prevent it from fragmenting into particles until it 
is well outside the nucleus. This feature immediately eliminates the possibility of a cascade 
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in the nucleus from the rescattering of the secondary products. If one makes the further 
assumptions that an excited nucleon interacts with the same cross section as an unexcited 
nucleon and that the successive collisions of the excited nucleon do not affect the excited 
state or its eventual fragmentation products [13], this leads to the conclusion (c. 1977) that 
the elementary process for particle production in nuclear collisions is the excited nucleon, 
and to the prediction that the multiplicity in nuclear interactions should be proportional to 
the total number of projectile and target participants, rather than to the total number of 
collisions, or R A = (N ch ) hA / {N ch ) hp = (N part ) hA / (N part ) hp = (1 + u)/2, as observed. This 
is called the Wounded Nucleon Model (WNM) [14] and, in the common usage, Wounded 
Nucleons (WN) are called participants. In a later model from the early 1980’s, the Additive 
Quark Model, AQM |T5], constituent-quark participants were introduced; but the AQM is 
actually a model of particle production by color-strings in which only one color-string can be 
attached to a constituent-quark participant, effectively a projectile quark participant model. 


5.2 Extreme Independent Models 


The models mentioned above are examples of Extreme Independent Models in which the 
effect of the nuclear geometry of the interaction can be calculated independently of the dy¬ 
namics of particle production which can be taken directly from experimental measurements. 
The nuclear geometry is represented by the relative probability, w n per A+B interaction 
for a given number n of fundamental elements, in the present case, number of collisions 
(Aeon), number of nucleon participants (wounded nucleon model-WNM [14]), number of 
constituent-quark participants (N qp ), number of color strings (AQM). The dynamics of par¬ 
ticle production, the AA or E?- distribution of the fundamental element, is taken from the 
measured p+p data in the same detector: e.g. the measured A c h distribution for a p+p colli¬ 
sion represents: 1 collision; 2 participants (WNM); a predictable convolution of constituent- 
quark-participants (NQP), or projectile-quark-participants (AQM). Glauber calculations of 
the nuclear geometry (w n ) together with the p+p measurement provide a prediction for the 
A+B measurement in the same detector as the result of particle production by multiple 
independent fundamental elements. 

I became acquainted with these models in my first talk at a Quark Matter conference 
(QM1984) where I presented measurements of transverse energy distributions from p+p and 
a+a interactions at ^/s~^ N =31 GeV at the CERN-ISR (Fig. 11a [18. 19]). The transverse 
energy, E r , is a multiparticle variable defined as the sum 


Et = ' S ^ / E i sin#* dEiT(rj)/dr) = sin 9{rj) dE(rj) / dp, (2) 

i 


where 0 is the polar angle, rj = — In tan 6/2 is the pseudorapidity, E r is by convention taken 
as the kinetic energy for baryons, the kinetic energy + 2 for antibaryons, and the total 
energy for all other particles, and the sum is taken over all particles emitted into a fixed 
solid angle for each event. 

The transverse energy was introduced by high energy physicists mm as an improved 
method to detect and study the jets from hard-scattering compared to high pr single particle 
spectra by which hard-scattering was discovered in p+p collisions and used as a hard-probe 
in Au+Au collisions at RHIC. However, it didn’t work as expected: E t distributions, like 
A c h distributions, are dominated by soft particles near (pr) (e.g. see Ref. [8] for details). 
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Figure 11: (a) E T distributions in p+p, a+a HE] at ^s NN =31 GeV, with AQM and WNM calculations [2T| . 
(b),(c) Ej- distributions with breaks indicating jets: (b) p+p y/s =62.3 GeV [25]; (c) da/dEr (nb/GeV) vs. 
E t for p+p yfs =540 GeV [23] . 


Nevertheless, it was claimed at the conference Em, in comments to my talk, that the devi¬ 
ation from the WNM in Fig. 0 was due to jets, but in both proceedings mm it was 
demonstrated that [20J “there is no ... sign of jets. This indicates that soft processes are 
still dominant, and that we are still legimately testing the WNM at these high values of E T .” 
As shown in Fig. Ink, the the AQM EH, rather than the WNM, followed the data. Jets do 
appear in E r distributions as a break <10 -5 down in cross section (Figs. 11 o, c). 

5.3 E t and N c h distributions cut on centrality 


At RHIC, following the style of the CERN SpS rather than the BNL-AGS fixed target 
heavy ion program, E T and N c h distributions were not generally shown. The measurements 
were presented cut in centrality in the form l^dNp^^ / dp /((N part ) /2) vs. (Ap art ) (Fig. [l2|), 
which would be a constant equal to ( dN c h pp /dr /) if the WNM worked. The measurements 
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Figure 12: (a) PHENIX, Au+Au, ^sj ~ =130 GeV [H]; (b) ALICE, Pb+Pb, =2.76 TeV jH]; (c) 

PHENIX preliminary Au+Au, ~ ~ =7.7 GeV compared to the data at larger ~ ~ [27 . 


clearly deviate from the WNM (Fig. |I2^ ) [21]; so the PHENIX collaboration, inspired by 
the preceding article in the journal [25] . fit their data to the two-component model: 

dN ch AA /d V = (dN ch pp /d V ) [(1 - x) (A^ part ) /2 + x (Won)] (3) 


where the N co \\ term implied a hard-scattering component for E T and 7V c h, known to be ab¬ 
sent in p+p collisions^] (recall Fig. 11). A decade later, the first measurement from Pb+Pb 

2 It was noted in Ref. [24 that hard-scattering was not a unique interpretation. The shape of the centrality 
dependences of d+T AA /dp parameterized as N part a were the same for Pb+Pb at +s NN =17.6 GeV at the 
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collisions with =2.76 TeV at the LHC Fig. 12 j [26] , showed exactly the same shape vs. 

the RHIC Au+Au data at ^/$ NN =200 GeV, although (N co n) is a factor of 1.6 larger 


AT, 


part 

and the hard-scattering cross section is more than a factor of 20 larger. This strongly argued 
against a hard-scattering component and for a nuclear geometrical effect, which was rein¬ 
forced by a PHENIX preliminary measurement in Au+Au at A+Cv =7.7 GeV (Fig. 12:) [27] 
which also showed the same shape for the evolution of (diV c h AA /dp /((Apart) /2) with A part 
as the sJ$ NN =200 and 2760 GeV measurements. It had previously been proposed that the 
number of constituent-quark participants provided the nuclear geometry that could explain 
the RHIC Au+Au data without the need to introduce a hard-scattering component [2BJ- 
However an asymmetric system is necessary in order to distinguish the NQP model from the 
AQM so the two models were applied to the RHIC d+Au data. 


5.4 The number of constituent-quark participants model (NQP) 

The massive constituent-quarks |29], which form mesons and nucleons (e.g. a proton =uud), 
are relevant for static properties and soft physics with p T <lA GeV/c. They are complex 
objects or quasiparticles Pi made of the massless partons (valence quarks, gluons and sea 
quarks) of DIS [31J such that the valence quarks acquire masses ~ 1/3 the nucleon mass with 
radii ~ 0.3 fm when bound in the nucleon. With smaller resolution one can see inside the 
bag to resolve the massless partons which can scatter at large angles according to QCD. At 
RHIC, hard-scattering starts to be visible as a power law above soft (exponential) particle 
production only for p T > 1.4 GeV/c at mid-rapidity (Fig[l3^), where Q 2 = 2 jjj, = 4 (GeV/c) 2 
which corresponds to a distance scale (resolution) <0.1 fmj^] 




Figure 13: (a) Invariant cross section of i r° vs. pr at mid-rapidity in p+p collisions at a/s =200 GeV [34] . 
The inset shows the transition from an exponental to a power-law in the range 1 < Pt < 2 GeV/c (b) 
PHENIX deconvolution of p+p E t distribution at y/s =200 GeV [1] 

CERN SpS and Au+Au at yjs NN =130 GeV, with a = 1.1 and a = 1.16 + 0.04, respectively. The LHC data 
10 years later [26] gave a = 1.19 ± 0.02 for Pb+Pb at . sJs NN =2760 GeV, again the same shape. 

3 Shuryak and collaborators recently made similar arguments about resolution in separating hard from soft 
processes although their mechanism for soft particle and QGP production was 2 color strings per wounded 
nucleon from Pomeron exchange [32] . 
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A standard Monte Carlo Glauber calculation is used to assemble the initial positions 
of all the nucleons. Then three quarks are distributed around the center of each nucleon 
according to the proton charge distribution p{r) oc e -VL 2 r/r rms , where r rms = 0.81 fm is the 
rms charge radius of the proton [33J. The q-q inelastic scattering cross section is adjusted to 
9.36 mb at s/s =200 GeV to give the correct p+p inelastic cross section (42 mb) and then 
used in the A+B calculations. 




Figure 14: PHENIX NQP calculations [F based on the Ex distribution of a constituent-quark participant 
from Fig. 13 ) for: (a) d+Au (also AQM), (b) Au+Au Ex distributions at s~ =200 GeV. 


Fig. [T3} ) shows the deconvolution of the p+p Ex distribution to the sum of 2-6 constituent- 
quark participants from which the Ex distribution of a constituent-quark is determined and 
applied to d+Au (Fig. [Tdji) and Au+Au (Fig. [Td})) reactions in the same detector. 

The NQP calculations closely follow the measured d+Au and Au+Au E T distributions in 
shape and magnitude over a range of more than 1000 in cross section. A complete calculation 
was also done for the AQM which fails to describe the d+Au data (Fig.[l4^). The conclusion 
is that the number of constituent-quark participants determines the N c ^ and Ex distributions 
and that the AQM calculation which describes the a-a data at y/s NN =31 GeV (Fig. Ill) 
is equivalent to the NQP in the symmetric system. 

The NQP model was also applied to the centrality-cut PHENIX data by making a plot 
of dEx/dr/ for a given centrality bin as a function of the number of constituent-quark par¬ 
ticipants N qp in the bin for Au+Au collisions at = 62.4, 130 and 200 GeV (Fig. 15 1 ). 



d) 

0 


o 3- 


O O 


PH ENIX 

preliminary 


O 2.76 TeV Pb+Pb (ALICE) 

• 200 GeV Au+Au 

■ 130 GeV Au+Au 

V 62 GeV Au+Au 

• 39 GeV Au+Au 
+ 27 GeV Au+Au 
A 19.6 GeV Au+Au 

• 7.7 GeV Au+Au 




IMMM * * 


♦ 


• • • 


b) 
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The data for each are well described by a straight line and all are consistent with a 

zero intercept. This means that dE T /dr] is strictly proportional to N qp (Fig [l5j:i) so that 
dET/drj/(N qp /2) vs. N qp is a constant for ^/s NN > 39 GeV (Fig [153) even up to the LHC 
^/s NN =2.76 TeV. This brought up a very interesting question, with a very important and 
fundamental answer. 

Most experiments at RHIC, starting with PHENIX (Fig. 12) [24] had successfully fit their 
measurements of dE T /dr] or dN^Jdr/ as a function of centrality (represented by iV part ) to 
the two-component model (Eq. [3j). Also, both the ATLAS [25] and ALICE [2,6] experiments 
at the LHC computed the ansatz, [(1 — x) (A r part ) /2 + x (A r co ii)], in event-by-event Glauber 
Monte Carlo calculations which fit their forward Er measurements used to define centrality 
in Pb+Pb collisions. ALICE realized that the combination of the two components Ap art 
and A coll in the ansatz represented the number of emitting sources of particles, which they 
named “ancestors”. Since the number of constituent-quarks N qp also represents the number 
of emitting sources in a simple linear relationship, Bill Zajc of PHENIX suggested that “the 
success of the two component model is not because there are some contributions proportional 
to A/p art and some going as A co n, but because a particular linear combination of A part and 
Aeon turns out to be an empirical proxy for the number of contsituent-quarks.” We had a 
nice table of (A pa rt), {N qp ), (N co \\) as a function of centrality in Au+Au collisions at ^ys~^ r 
=200 GeV, so it did not take very long for me to verify the striking result that indeed 
it was true: with x = 0.08, the ratio (N qp ) /[(1 — x) (A part ) /2 + x (A co u)] = 3.88 on the 
average and varies by less than 1% over the entire centrality range in 5% bins except for the 
most peripheral bin where it is 5% low (Table [l]). This result clearly demonstrates that the 
ansatz works because the particular linear combination of A part and A co n turns out to be an 
empirical proxy for N qp and not because the iY co u term implies hard-scattering. 




Centrality 

(Npart) 

(Nqp) 

(iVcon) 

ansatz 

(7V q p)/ansatz 

0 - 5 % 

350.9 ±4.7 

956.6 ± 16.2 

1064.1 ± 110.0 

246.5 

3.88 

5 - 10 % 

297.0 ±6.6 

789.8 ± 15.3 

838.0 ± 87.2 

203.7 

3.88 

10-15% 

251.0 ±7.3 

654.2 ± 14.5 

661.1 ±68.5 

168.3 

3.89 

15-20% 

211.0 ±7.3 

540.2 ± 12.3 

519.1 ±53.7 

138.6 

3.90 

20-25% 

176.3 ±7.0 

443.3 ± 10.4 

402.6 ± 39.5 

113.3 

3.91 

25-30% 

146.8 ±7.1 

362.8 ± 12.2 

311.9 ±31.8 

92.5 

3.92 

30-35% 

120.9 ±7.0 

293.3 ±11.0 

237.8 ± 24.2 

74.6 

3.93 

35-40% 

98.3 ±6.8 

233.5 ± 9.2 

177.3 ± 18.3 

59.4 

3.93 

40-45% 

78.7 ±6.1 

182.7 ± 6.8 

129.6 ± 12.6 

46.6 

3.92 

45-50% 

61.9 ±5.2 

140.5 ± 5.3 

92.7 ±9.0 

35.9 

3.91 

50-55% 

47.6 ± 4.9 

105.7 ± 5.5 

64.4 ±8.1 

27.0 

3.91 

55-60% 

35.6 ±5.1 

77.3 ± 6.8 

43.7 ± 7.6 

19.9 

3.89 

60-65% 

26.1 ±4.7 

55.5 ±7.1 

29.0 ±6.5 

14.3 

3.87 

65-70% 

18.7 ±4.0 

39.0 ±6.7 

18.8 ± 5.3 

10.1 

3.86 

70-75% 

13.1 ±3.2 

27.0 ± 4.9 

12.0 ± 3.6 

7.0 

3.86 

75-80% 

9.4 ±2.1 

19.0 ±3.2 

7.9 ± 2.2 

5.0 

3.83 

80-92% 

5.4 ±1.2 

10.3 ±1.5 

4.0 ±1.0 

2.8 

3.67 

P+P 

2 

2.99 ± 0.05 

1 

1 

2.99 




Table 1: Verification that the ansatz, [(1 — x) (iV part ) /2 + x (iV co n)], from Eq. [5] with x = 0.08, is a proxy 
for Ngp. The errors quoted on (N part ), (Ag p ), (N co n) are correlated and largely cancel in the {N qp }/ ansatz 
ratio. For x = 0.09 the average (N qp ) /ansatz=3.81, the maximum variation is less than 1.6%, but 4% low 
in the most peripheral bin. 
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The fact that the (N qp ) /ansatz ratio drops from an average of 3.88 for Au+Au collisions 
to 2.99 for p+p collisions is also interesting. This is consistent with the PHOBOS Ea result 
that a fit of Eq. rh to dN c + A /dr /) with x = 0.09 leaving (dN c h pp / drj) as a free parameter 
gives the result [dN c ] 1 pp / drj) = 2.70 which is above the measured inelastic value of 2.29. 
The lower value of (N qp ) /ansatz for p+p would then give a value of 2.70 x 2.99/3.88 = 2.08 
(2.12 for x = 0.09) for {dN c \ pp /dr/) , much closer to the measured value. In that same paper, 
PHOBOS also noted that their data were consistent with a constant value of x from ^/$ NN 
=19.6 to 200 GeV (more recently extended to ^$ NN =2.76 TeV [BE]) which indicated that 
the fraction of hard-processes contributing to multiplicity did not increase over a huge range 
of ■s/s even though the hard-scattering cross section greatly increased over this same range. 


5.5 Constituent-quark participants resolve several outstanding puz¬ 
zles 

PHOBOS also made some very nice measurements of the charged particle multiplicity over 
the full rapidity range, not just mid-rapidity. The total charged multiplicity iV c h||, ) |< 5.4 
was measured and 2V c h| |tj|< 5 . 4 / (A/ )art /2) plotted as a function of centrality for ^/$ NN =19.6, 
62.4,130, and 200 GeV (Fig. [Toj -i) [39]. At first glance the data appear to follow the WNM 
because the multiplicity/per nucleon pair A c h|| ?7 |< 5 . 4 / (N paxt /2) appears to be constant in 
Au+Au collisions. However, the true believers, e.g. Ref. m, claim that the WNM does not 
work because the value in Au+Au collisions is much larger than the p+p value shown; but 
“still the proportionality of these multiplicities to the number of participants holds” [TO]. 
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Figure 16: a) Total charged multiplicity per nucleon pair N^/ (Af part /2) vs. centrality, Af part , for the ^s NN 
indicated. Open points are the measured iV c h|| 77 |<5.4/((iV part /2); solid points are extrapolated to \r]\ < 2/beam- 
b) Total charged multiplicity per nucleon pair in p+p and A+A collisions as a function of c.m. energy yfs 
compared to e + + e~ collisions mi- 
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In fact, the difference between the p+p and Au+Au values may be related to another 
interesting observation by PHOBOS m that the “leading particle effect” in p+p collisions, as 
discovered by Zichichi and collaborators [42]—in which the total multiplicity at c.m. energy 
WsW is equal to that in e + e _ collisions at = 


(the “effective energy”) (Fig. 16 d) 


because the leading protons carry away half the p+p c.m. energy—is absent in A+A collisions 
where the leading protons can reinteract. This observation seems to contradict the WNM, 
in which the key assumption is that what counts is whether or not a nucleon was struck, not 
how many times it was struck. 

Both these effects can be reconciled by constituent-quark participants. 

In the NQP model (Table[l|, the {N qp /N part ) is 1.5 for a p+p collision but rises to 2.27- 
2.73 (a factor of 1.51-1.82) for the more central (0-50%, (7V part ) > 60) Au+Au collisions 
plotted in Fig. 161. This would correspond to an increase in A%/ {N paxt /2) by a factor of 
~ 1.5 as observed. It also might explain the slight rise of the open points with increasing 
Afpart- Similarly, the increase in “effective energy” for particle production shown by the 
increase in A%/ (A part /2) from p+p to Au+Au collisions is due to an increase in the number 
of (constituent-quark) participants, not because of additional collisions of a given nucleon- 
participant. Furthermore, the factor 1.5 decrease in A%/ (A part /2) from Au+Au to p+p 
corresponds to a reduction in d/s for the observed N c ^/ (N part /2) from 200 to 100 GeV 
p+p collisions on Fig. 16 + the same factor of 2 discussed in the original measurement [TlJ. 
Thus, the NQP model rather than the WNM preserves the assumption in these “extreme- 
independent” participant models that successive collisions of a participant do not increase 
its particle emission while explaining these two interesting observations. 

Another argument against the -A part , A+.,h ansatz representing actual hard-collisions 
rather than simply being a proxy for constituent-quark participants concerns the measure¬ 
ment of elliptic flow in central U+U collisions 


6 Collective Flow 

For many years, since the days of the Bevalac & collective flow 03 has been observed in 
A+A collisions over the full range of energies studied, from incident kinetic energy of 100A 
MeV to c.m. energy of d/s// = 2.76 TeV, and thought to be a distinguishing feature of 
A+A collisions compared to either p+p or p+A collisions. Collective flow, or simply flow, 
is a collective effect which can not be obtained from a superposition of independent N+N 
collisions. I first present a short review (details can be found in previous ISSP proceedings [3 
m and then move on to the newer results. 

Immediately after an A+A collision, the overlap region defined by the nuclear geometry 
is almond shaped (see Fig |I7^ ) with the shortest axis along the impact parameter vector. 
The different pressure gradients along the short and long axes of the ellipse break the 0 
symmetry of the problem and create an azimuthal angular dependence of the semi-inclusive 
single particle spectrum with respect to the reaction plane, 0 — <Fr, which is represented by 
an expansion in harmonics [48], where the angle of the reaction plane is defined to be 
along the impact parameter vector, the x axis in Fig. 

Ed 3 N _ d 3 N _ cl 3 N 
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Figure 17: (a) Almond shaped overlap zone generated just after an A+A collision where the incident nuclei 
are moving along the ±z axis. The reaction plane by definition contains the impact parameter vector (along 
the x axis) [40]. (b) v 2 for charged particles integrated over pr at y/s NN =2.76 TeV for 20-30% centrality 
compared to the measurements at lower ^/s NN at the same centrality [471 . 


The Fourier coefficient v 2 , called elliptic flow, is predominant at mid-rapidity. The evolution 
of v 2 with y/s^ (Fig. 17 d) [U] is the result of competing processes. At very low 
corresponding to values of ~ 100A MeV |49j the main effect among many others is from 
nuclei bouncing off each other and breaking to fragments, which is sensitive to the equation 
of state of the nuclei—soft, like sponges, hard like billiard balls? The negative v 2 at larger 
v /s^ is produced by the effective “squeeze-out” (in the y direction) of the produced particles 
by slow moving minimally Lorentz-contracted spectators (as in Fig. [TTf t) which block the 
particles emitted in the reaction plane. With increasing the spectators move faster 

and become more contracted so the blocking stops. The increase of v 2 with ^/s NN is generally 
described by hydrodynamics in the QGP region, but is also described by hadron transport 
theories for • v /s JVJV <10 GeV [HD] , 

Flow measurements contributed two of the most important results about the properties 
of the QGP: i) the scaling of v 2 of identified particles at mid-rapidity with the number of 
constituent-quarks n q in the particle— v 2 /n q scales with the transverse kinetic energy per 
constituent-quark, KE T /n q , because particles have not formed at the time flow develops; 
ii) the persistence of flow for Pt > 1 GeV/c which implied that the viscosity is small [fLlj , 
perhaps as small as a quantum viscosity bound from string theory [52], p/s = 1/(47 t) where 
p is the shear viscosity and s the entropy density per unit volume. This led to the description 
of the “sQGP” produced at RHIC as “the perfect fluid”. 

New insight came in 2013, when measurements in p+Pb at LHC and d+Au at RHIC 
observed what looked very much like collective flow in these systems that were believed 
to be too small to support collective effects. This was the reason for the He 3 +Au run at 
RHIC in 2014, to see whether triangular flow, v 3 , would be more prominent with a 3 nucleon 
projectile. The improvement of the p+Au and ci+Au measurements this year to identified 
pions and protons strengthened the case that the observed v 2 in these small systems is really 
hydrodynamic collective flow. 

Figure 18 1 [53J shows the two-particle correlation function in rf+Au (Eq. [5]) fit with 
terms from C\ to c 4 where the solid line is the fit and only the cp (dashes) and c 2 (dots) make 
significant contribution. 
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Figure 18: (a) PHENIX [53] two-particle azimuthal correlation function (Eq. [5]) for 1.0 < PT,trig < 3.0 
GeV/c, |?7trig| < 0.35, in central (0-5%) d +Au collisions at RHIC. ^(pt) for identified p ± by the standard 
reaction plane method (Eq. [4| for (b) central d +Au collisions at RHIC (y/s NN =200 GeV) and (c) central 
p+Pb collisions at LHC (Ws NN =5.02 TeV). 
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The ci comes into play because the trigger particle t is a charged track at mid-rapidity while 
the associated particle a is a count in an MPC tower (Sr) x S(p ~ 0.12 x 0.18) from 7 T° or 
7) meson decay photons at —3.7 < r/ tiOW er < —3.1. Also, there is no evidence of a di-jet 
contribution because the large pseudorapidity gap between t and a is beyond that of a di-jet. 
Thus, the long-range correlation in Fig. [18^ which is not seen in p+p comparison data but has 
the same properties as collective flow in Au+Au collisions is consistent with hydrodynamic 
collective flow in cJ+Au. Perhaps more convincing evidence for hydrodynamic flow is given 


in Fig. 18 3 ,c where both at RHIC in d+Au (b) and LHC in p+Pb (c), the characteristic 7 r, 
p mass splitting for v 2 (pt ) seen in Au+Au is observed [ 53 i. The splitting occurs because, 
for a given transverse collective expansion velocity (3, protons have a larger pr = 7 /3m than 
pions. 

6.1 V 2 in U+U collisions and constituent-quark participants 

Because Uranium nuclei are prolate spheroids, there is the interesting possibility of large v 2 
in body-to-body central collisions which have a significant eccentricity and almond shape 

ansatz (Eq. [3]) would describe 


(Fig. 19 a). Based on the assumption that the iV part , N co \\ 


the dN c h/dri distribution in U+U collisions, it was predicted that for the highest dN^-Jdrj 
(the most central collisions) the tip-to-tip configuration with much larger N co u and small 
eccentricity (small v 2 ) would overtake the body-to-body configuration with large eccentricity 
corresponding to large v 2 . 

This led to two predictions: i) the tip-to-tip configuration would be selected by the most 
central collisions [Mj; ii) these most central collisons would see a sharp decrease in v 2 with 
increasing dN^/dr) [ 55 ], 56| called a cusp. This sharp decrease—represented by the bent 
line on the topmost U+U data (filled circles) in Fig. [T 9 } ) (called a knee in Ref. [43])—is 
not observed. As discussed previously, this is because the N co \\ term is not relevant for 
d,N c \Jdr) distributions, which also argues against the method proposed in Ref. [ 54 J to select 
the tip-to-tip configuration. 
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Figure 19: (a) Body-to-body and tip-to-tip configurations in U+U collisions with zero impact parameter. 
The different relation of 7V part to 7V co n is sketched next to each configuration. (Modified drawing from 
Ref. [43]). (b) STAR measurements of V 2 in Au+Au and U+U at ^/s NN « 200 GeV as a function of dN^/dr] 
with upper percentiles of centrality for U+U indicated by vertical dashed lines [43] . 


7 RHIC Beam Energy Scan (BES)—in search of the 
critical point 

In addition to discovering the QGP and measuring its properties, another objective of the 
RHIC physics program is to measure the phase diagram of nuclear matter and to determine 
the equation of state in the various phases and the characteristics of the phase transitions. 
Two of the many proposed phase diagrams of nuclear matter (e.g. see Ref. EH) are shown 
in Fig. [20] together with the idealized trajectories of the evolution of the medium for Au+Au 
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Figure 20: Proposed phase diagrams for nuclear matter: Temperature, T, vs Baryon Chemical Potential, 
He- a)STAR’s idea in 2013 [SB]; b) STAR’s more cautious idea in 2014 [59]. 


collisions at the y/s~^ proposed for the Beam Energy Scan at RHIC to search for a QCD 
critical point. The bursts represent the hottest and densest stage of the medium when 
thermal equilibrium is reached shortly after the collision. The axes are the temperature T 
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vs. the baryon chemical potential ps■ The temperature for the transition from the Quark 
Gluon Plasma (QGP) to a hadron gas is taken as 170 MeV for ps = 0 and the phase 
boundary is predicted to be a smooth crossover down to a critical point below which the 
phase boundary becomes a first order phase transition. 

In an equilibrated thermal medium, particles should follow a Boltzmann distribution in 
the local rest frame [60] 


d 2 a 

dpLPrdpr 


d 2 a 


oc 


dpLrriTdmT e*T G/ T ± 1 


~ e 


-(E-u,)/T 


( 6 ) 


where mp = \JpT + 11,2 anf l P is a chemical potential. In fact, the ratios of particle abun¬ 
dances (which are dominated by low px particles) for central Au+Au collisions at RHIC, 
even for strange and multi-strange particles, are well described [61] by fits to a thermal 
distribution, 


_ ^ p -(E-n)/T P = e_ = _(2 hb)/t ( 7 \ 

dpLPrdpr P 

with similar expressions for strange particles, ps (and p$) are chemical potentials associated 
with each conserved quantity: baryon number, pb, (and strangeness, ps). Thus it is simple 
and instructive to estimate the p/p ratio from Fig. [20[ :t, near the arrow, where I read T = 160 
MeV, ps = 300 MeV, ~ 30 GeV, which gives p/p ~ 0.02. Since the p/p ratio vs 

will be an important issue later, it is not a good idea to get this important information from 
a sketch in a proposal but from measurements and the best analysis (Fig. J2l|. The results 




Figure 21: (a) STAR measurements of p/p vs. \f* NN 157]: b) Best accepted analysis of T and pb vs 

VTG7 [S3- 


are: 


i) the correct p/p ratio at y/s//// = 30 GeV is 


0.12 from Fig. 21 1 ; ii) this ratio also 
corresponds to the correct pb ~ 170 MeV at ^/$ NN = 30 GeV from Fig. [2TJd. The lesson is: 
if it looks more like art than like science, be skeptical and look in refereed journals for the 
correct numbers. 
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7.1 A press release during ISSP 2011 

On June 23, 2011, shortly before I was to present my 2011 lectures, a press release from LBL 
arrived claiming that “By comparing theory with data from STAR, Berkeley Lab scientists 
and their colleagues map phase changes in the QGP” [64]. Since I was going to criticize 
in my lectures what I considered to be a particularly egregious case of “physics by press 
release” in the year 2000 by CERN (see Ref. !),i felt that I was obliged to quickly absorb 
and present in my talk the physics behind this latest press release, hopefully a “Highlight 
from RHIC”. 

The subject is “Fluctuations of conserved quantities”, in this case the net baryon distri¬ 
bution taken as p — p. Since there can be no fluctuations of conserved quantities such as 
net charge or net baryon number in the full phase space, one has go to “locally conserved 
quantities” [HHJ in small rapidity intervals to detect a small fraction of the protons and anti¬ 
protons which then fluctuates, i.e. varies from event to event. The argument is that, e.g. 
the fluctuation of one charged particle in or out of the considered interval produces a larger 
mean square fluctuation of the net electric charge if the system is in the hadron gas phase 
with integral charges than for the QGP phase with fractional charges. 

However, while there are excellent statistical mechanical arguments about the utility 
of fluctuations of conserved quantities such as net baryon number as a probe of a critical 
point [66], there were, in 2011, no adequate treatments of the mathematical statistics of the 
experimental measurements. There are also additional problems such as short-range rapidity 
correlations in A+A collisions between like-particles induced by Fermi or Bose quantum 
statistics that must be reckoned with (e.g. see Refs. mm)- 

Theoretical analyses tend to be made by a Taylor expansion of the free energy F = —T In Z 
around the critical temperature T c where Z is the partition function, or sum over states, 
which is of the form 

Z oc e -(£-EiftQi)AP 

and fa are chemical potentials associated with conserved charges Q - t [66]. The terms of the 
Taylor expansion, which are obtained by differentiation, are called susceptibilities, denoted 
X- The only connection of this method to mathematical statistics is that the Cumulant gen¬ 
erating function in mathematical statistics for a random variable x is also a Taylor expansion 
of the In of an exponential: 



9x(t) = In ( e tx ) = V kJ— 

'' n\ 


ftm, — 


d m g x (t ) 


n— 1 


dt Ti 


(9) 


t =0 


Thus, the susceptibilities are Cumulants in mathematical statistics terms, where, in general, 
the Cumulant n rn represents the m th central moment, = ((x — i-i) rn ), with all m-fold 
combinations of the lower order moments subtracted, where // = ( x }. For instance, hy = 
((x — h) 2 ) = cr 2 , K 3 = ((x — h) 3 )i k 4 = (i x ~ h) 4 ) — 3 ^ 2 , «5 = ((x ~ h) 5 ) ~ Wk,sk, 2 - Two 
so-called normalized or standardized Cumulants are common in this field, the skewness, 
S = K 3 / cr 3 and the kurtosis, k = fiq/cr 4 = {(x — /i) 4 ) /a 4 — 3. 

A sample [69] of STAR measurements of the distribution of net-protons in Au+Au col¬ 
lisions in the small interval 0.4 < p T < 0.8 GeV/c, y\ < 0.5 for different QTvT is shown in 
Fig. 22 a. The moments in the form na 2 = k.i/k 2 are shown from a previous STAR publica¬ 


tion ua in Fig. |22p while a plot, alleged to be of this same data, presented in the Lattice 
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Figure 22: a) (top-left) STAR |69] distribution of event-by-event p — p at 3 values of y/ snn ; b) (top-right) 
STAR published [70] measurements of kg 2 ] c) (bottom-left) Measurements from (b) as shown in Ref. [7T] 
compared to the predicted ratio of susceptibilities (open crosses); d) (bottom-right) compilation [6l| of STAR 
measurements of kg 2 for p — p. 


QCD theory publication that generated the press-release, is shown in Fig. 22 
plot of the kg 2 from the data of Fig. 22 3, combined with the results from Fig. 


in Fig. 221 [69] . There are many interesting issues to be gleaned from Fig. p22 


[71] ; and 


22 o, is shown 


The data point at 20 GeV in Fig. 22: is not the published one from (b), as stated in the 
caption im but the one from (d), which is different and with a much larger error. This, in 
my opinion, makes the data point look better compared to the predicted discontinuous value 
of ncr 2 = —1.5 for the critical point at 20 GeV (open crosses) in contrast to the predictions 
of 1.0 for both 62.4 and 200 GeV. The published measurements in (b) together with the 
newer measurements in (d) are all consistent with kct 2 = 1; but clearly indicate the need for 
a better measurement at ysjvw = 20 GeV. Apart from these issues, the main problem of 
comparing Lattice QCD “data” to experimental measurements is that it is like comparing 
peaches to a fish, since the prediction is the result of derivatives of the log of the calculated 
partition function of an idealized system, which may have little bearing on what is measured 
using finite sized nuclei in an experiment with severe kinematic cuts. Maybe this is too harsh 
a judgement; but since this is the first such comparison (hence the press release), perhaps 
the situation will improve in the future. If a future measurement would show a significant 
huge discontinuity of na 2 similar to the theoretical prediction at yG/vw = 20 GeV, then even 
I would admit that such a discovery would deserve a press release, maybe more! 
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7.1.1 If you know the distribution, you know all the moments and cumulants 


When I first saw the measured distributions in Fig. [22 p in 2011, my immediate reaction was 
that STAR should fit them to Negative Binomial distributions (NBD) so that they would 
know all the Cumulants. However, I subsequently realized that my favorite 3 distributions for 
integer random variables, namely, Poisson, Binomial, and Negative Binomial, are all defined 
only for positive integers (e.g. see Ref. [E] for details), while the number of net-protons on 
an event can be negative as well as positive, especially at higher c.m. energies. Thanks 
to Gary Westfall of STAR, in a paper presented at the Erice School of Nuclear Physics 
in 2012 [72], who found out that these three distributions fall into the class of “integer 
valued Levy processes ra” for which the Cumulants Hj for the distribution P(n — m ) of the 
difference of samples from two such distributions, P + (n ) and P~(m ), with Cumulants kJ 
and kJ, respectively, are HUES]: 


= K j + ( _ 1 ) J/ S- » ( 10 ) 

so long as the distributions are not 100% correlated. This result is the same as if the 
distributions P + (n ) and P~{m ) were statistically independent. The first four Cumulants of 
the Poisson, Binomial and Negative Binomial distributions are given in Table [2j 


Table 2: Cumulants for Poisson, Binomial and Negative Binomial Distributions 


Cumulant 

Poisson 

Binomial 

Negative Binomial 

K\— 11 

H 

np 

p 

«2 = /A = <J 2 

/-' 

p{i-p) 

p( 1 + p/k) 

«3 = Hz 

H 

a 2 ( 1 - 2 p) 

cr 2 (l + 2 p/k) 

K 4 = /i4 — 3 ac| 

H 

cr 2 (l — 6p + 6 p 2 ) 

cr 2 ( 1 + Qp/k + 6 p 2 /k 2 ) 

S = k 3 /ct 3 

1/Vh 

(1 - 2 p)/a 

(1 + 2 p/k)/a 

K = K^/k 2 

1/p 

(1 — Qp + 6 p 2 )/a 2 

(1 + Qp/k + Qp 2 /k 2 )/a 2 

Sa = k 3 /k 2 

1 

(1 - 2 p) 

(1 + 2 p/P) 

na 2 = K 4 /k ,2 

1 

(1 — 6p + 6 p 2 ) 

(1 + Qp/k + Qp 2 /k 2 ) 


7.2 The latest measurements have appeared without a press re¬ 
lease. 

In the intervening period since 2011, the STAR collaboration has improved the preliminary 
measurements to publications and has improved the analysis by comparing to both Poisson 
and Negative Binomial distributions. Figure 23 a, [75] shows the STAR measurements of 


Cumulants of the net charge {N + — N~) distributions from the “number of positive (A r+ ) 
and negative ( N ~) charged particles within \rj\ < 0.5 and 0.2 < pr < 2.0 GeV/c on each event 
(after removing protons and antiprotons with pr < 400 MeV/c) [ZH]”- The corresponding 
Poisson and NBD Cumulants were calculated from the measured mean, p, and variance, 
cr 2 , of the N + and N~ distributions, respectively, and then calculated using Eq. 10, In 


contrast to Fig. 22 


no non-monotonic behavior with s/s^n is observed (or claimed) and 
the measurements of Sa and na 1 are all above the Poisson baseline. The Sa measurements 


clearly favor the NBD. 
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Figure 23: a) yj s^n dependence of combinations of Cumulants in Au+Au (and p+p) from STAR: a) (left) 
net-charge Cumulants m > where M is used to represent the mean, p. b) (right) Cumulants of the N. p 
distributions [76] , where the error bars are statistical and the caps systematic errors. 


N p 


The situation is quite different for the net-proton (N p — N p ) Cumulants (Fig. |'23 
measured within \y\ < 0.5 over the range 0.4 < p? < 0.8 GeV/c which covers roughly half 
the pt spectrum. Here the measurements of Sa and kg 2 are all below the Poisson baseline, 
denoted Skellam, which is the distribution of the difference between two Poissons and reflects 
“a system of totally uncorrelated, statistically random particle production” [76]. From Eq. 10 
and Table [2] for a Poisson one can see that, for a Skellam, Sa = K 3+2 = (l- f p ~ /%)/ (l J p + /%) 
which increases with decreasing ^/s NN because the p vanish (/ip = (N p ) <C p p ) so that the 
shape of the net distribution becomes dominated by the protons. This is easier to see in a 
plot of kg 2 vs pb with ^/$ NN indicated (Fig. 24 r) [59] which shows clearly that kg 2 starts 


dropping for yjs~ < 39 GeV where the p/p ratio drops below ~ 0.3 (recall Fig. 21). 
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Figure 24: (a) (left) kg 1 vs. y/s NN and ps [5^; b) (right) ( dN p /dy — dNp/dy) for top 5% centrality at 
AGS (Au+Au, y/s NN =4.9 GeV), SPS (Pb+Pb, y/s NN =17.2 GeV) and RHIC(Au+Au, y 's NN =200 GeV), 
beam rapidity in c.m. system y v = 1.6, 2.9, 5.4 [77] . 
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The errors are still too large to determine whether or not kg 2 stays constant below ^fs 


=19.6 GeV but are sufficient to clearly rule out the value of kg 2 = —1.5 at ^J~s 


NN 

20 GeV 

predicted in Fig. |22| [[71] which created the fuss in 2011. It is also important to point out that 
in addition to the vanishing of the anti-protons, the physics of the protons at mid-rapidity 
changes dramatically—the protons are no longer produced particles, which would conserve 
N p — Np, but are the participants and fragments from the colliding nuclei which move to mid¬ 
rapidity and eventially stop as is reduced from 200, to 17.2 to 4.9 GeV (Fig. 24 :>) (73. 
All these results indicate that the search for a QCD critical point at RHIC in the Beam 
Energy Scan (BES) in 2018-19 may not be as straightforward as originally assumed. 


8 Jet quenching, RHIC’s main claim to fame 

The gold-plated signature for the QGP since 1986 [78] has been the suppression of J /+ 
because the color potential between c, c quarks would be screened (Debye screening) by all 
the free color charges in the medium so that the c, c would not be able to bind to form the 
J/T. In fact the PHENIX experiment at RHIC was specifically designed to detect the J/ + 
at mid-rapidity at rest or with very low pr (where the screening effect would be the largest) 
via the decay J / + —» e + + e _ . J / + suppression was reportedly observed several times 
at the CERN SpS fixed target heavy ion program starting with NA38 in O+U collisions 
in 1989 [79] but was plagued with many problems. The principal physics problem is that 
the J/ + does not follow the standard hard-scattering pointlike scaling in A+B collisions, 
g J J^ = (7b-B) Q-cr 7VJV with a = 1, but is suppressed in cold nuclear matter (CNM) inp+A and 
A+B collisions, with a ~ 0.91. Thus, the ultimate discovery by NA50 in Pb+Pb collisions 
at yfs NN =17.2 GeV [80] in 1996-1998 was called “anomalous suppression” because it was 
below the CNM cross section dependence which was itself well below the hard-scattering 
pointlike scaling. 

In 1998 at the QCD workshop in Paris ED, i found what I thought was a cleaner signal 
of the QGP when Rolf Baier asked me whether jets could be measured in Au+Au collisions 
because he had made studies in pQCD [82] of the energy loss of partons, produced by hard- 
scattering “with their color charge fully exposed”, in traversing a medium “with a large 
density of similarly exposed color charges”. The conclusion was that “Numerical estimates 
of the loss suggest that it may be significantly greater in hot matter than in cold. This makes 
the magnitude of the radiative energy loss a remarkable signal for QGP formation” [82] • In 
addition to being a probe of the QGP the fully exposed color charges allow the study of 
parton-scattering with Q 2 <C 1-5 (GeV/c) 2 in the medium where new collective QCD 
effects may possibly be observed. 

Because the expected energy in a typical jet cone R = yj (Ary) 2 + (A0) 2 in central Au+Au 
collisions at y/s NN =200 GeV would be ivR 2 x l/2n x dE^/dry = R 2 / 2 x dE^/dry ~ 350 
GeV for R = 1, where the kinematic limit is 100 GeV, I said (and wrote [81] ) that jets can 
not be reconstructed in Au+Au central collisions at RHIC—still correct after 16 years. On 
the other hand, hard-scattering was discovered in p + p collisions at the CERN-ISR in 1972 
with single particle and two-particle correlations, while jets had a long learning curve from 
1977-1982 with a notorious false claim (e.g. see Refs. M), so I said (and wrote [SI]) that 
we should use single and two-particle measurements—which we did and it WORKED! The 
present solution for jets in A+A collisions (LHC 2010 and RHIC c.2014) is to take smaller 
cones, with 56 GeV in R = 0.4, 32 GeV in R = 0.3, 14 GeV in R = 0.2 at RHIC. 
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8.1 Jet quenching at RHIC — Suppression of high px particles 


The discovery at RHIC [83] that 7r°’s produced at large transverse momenta are suppressed in 
central Au+Au collisions by a factor of ~ 5 compared to pointlike scaling from p+p collisions 
is arguably the major discovery in Relativistic Heavy Ion Physics. For 7T° (Fig. 25 1 ) [84] the 
hard-scattering in p-\-p collisions is indicated by the power law behavior p^ n for the invariant 
cross section, Ed 3 a/dp 3 , with n = 8.1 ± 0.1 for pr > 3 GeV/c. The Au+Au data at a given 
Pt can be characterized either as shifted lower in pr by dp' T from the pointlike scaled p+p 
data at p' T = Pt + 8p' T , or shifted down in magnitude, i.e. suppressed. In Fig. 25 3, the 


suppression of the many identified particles measured by PHENIX at RHIC is presented as 
the Nuclear Modification Factor, Raa(pt ), the ratio of the yield of e.g. 7r per central Au+Au 
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Figure 25: a) (left) Log-log plot of invariant yield of i r° at yj snn = 200 GeV as a function of transverse 
momentum p T in p+p collisions multiplied by (Taa) for Au+Au central (0-10%) collisions compared to the 
Au+Au measurement [84]. Vertical arrow is for Raa{pt ), horizontal arrow for Sp' T . b) (right) Raa{pt ) for 
all identified particles so far measured by PHENIX in Au+Au central collisions at y / W/v = 200 GeV. 


collision (upper 10%-ile of observed multiplicity) to the point like-scaled p+p cross section at 
the same pr, where (Taa) is the average overlap integral of the nuclear thickness functions: 

(1 /N A a) d 2 Nl A /dp T dy 


Raa(pt ) — 




AA 


d 2 Cp P 


/dp T dy 


( 11 ) 


The striking differences of Raa(pt) in central Au+Au collisions for the many particles 
measured by PHENIX (Fig. [25):)) illustrates the importance of particle identification for un¬ 
derstanding the physics of the medium produced at RHIC. Most notable are: the equal 
suppression of 7r° and r\ mesons by a constant factor of 5 (Raa = 0.2) for 4 < pr < 15 
GeV/c, with suggestion of an increase in Raa for pr > 15 GeV/c; the equality of sup¬ 
pression of direct-single e ± (from heavy quark (c, b ) decay) and 7T° at £+> 5 GeV/c; the 
non-suppression of direct-q for pr > 4 GeV/c; the exponential rise of Raa of direct-q for 
Pt < 2 GeV/c [HS), which is totally and dramatically different from all other particles and 
attributed to thermal photon production by many authors (e.g. see citations in Ref. [85j). 
For pt> 4 GeV/c, the hard-scattering region, the fact that all hadrons are suppressed, but 
direct-q are not suppressed, indicates that suppression is a medium effect on outgoing color- 
charged partons likely due to energy loss by coherent Landau-Pomeranchuk-Migdal radiation 
of gluons, predicted in pQCD [82] , which is sensitive to properties of the medium. 
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One nice advantage that hard-scattering and high p T suppression have as a QGP probe 
compared to J/T suppression is that although there is a CNM effect, it is an enhancment 
rather than a suppression; and as far as is known, the enhancement, historically called the 
Cronin effect [Eb], only occurs for baryons and not mesons at RHIC energies. Figure 26 r 



Figure 26: Measurements of Raa of identified particles 7r , K* , , 7r° , <p as a function pt and centrality 
at -y/Sjvjv = 200 GeV [57]: a) (left) Au+Au; b) (right) d+Au. 

shows Raa in Au+Au for protons and mesons in the range 0.5 < pr < 6.0 GeV/c, where, in 
central collisions (0-10%), all the mesons are suppressed for p T > 2 GeV/c while the protons 
are enhanced for 2 < p T < 4 GeV/c and then become suppressed at larger p T . The d+Au 
results in Fig. 26 d show no CNM effect for the mesons, Raa ~ 1 out to pr = 6 GeV/c; 


while the protons show a huge enhancement (Cronin effect) in all centralities except for the 
most peripheral (60-88%). At present, there is no explanation of the proton enhancement in 
either Au+Au or d+Au collisions, so 7T° and rj are the favored hard-probes. 

8.2 ^Pt/Pt’ the fractional shift in the p^ spectrum 

After more than a decade of using the ratio Raa, we are now paying more attention to Sp' T /p' T , 


the fractional shift of the p' T spectrum, as an indicator of energy loss in the QGP Fig. 27 


For a constant fractional energy loss, which is true at RHIC in the range 6 < pr < 12 GeV/c 


as shown in Fig. 25 1 where the p+p reference and Au+Au measurement are parallel on a 


log-log plot) there is a simple relationship between Raa, $ Pt/Pt an( l n , the power in the 
invariant pr spectra: 

Raa(p't ) = Raa(pt) = (1 — fi'jP t/ptY 2 • (12) 

Using Spr/pr is important for comparison to the LHC measurements where the power is 
n ~ 6 compared to n = 8.1 at RHIC, so that the same Raa does not mean the same Spr/pr- 
Strictly Spr/pr is not a measure of the parton energy loss in the QGP but is used as a 
proxy. Figure 27 1 shows that Spr/pr at pr = 7 GeV/c for RHIC and LHC both increase 


monotonically with centrality (A% irt ) but is a factor of 2 to 1.4 larger at LHC, depending 
on centrality, a likely indication of a hotter and/or denser medium. Figure [27)) attempts to 
determine whether Spr/pr is a universal function of the charged particle density, dN^/dp at 
both RHIC and LHC. The dependence is not quite universal. A fit of Spr/pr oc (dN^/dp) a 
gives a ~ 0.35 at LHC and 0.55 at RHIC, although the data at ^/s NN = 200 GeV and 
2.76 TeV do appear to merge for (dN^/dp) > 300. Hopefully, measurements of Spr/pr will 
eventually lead to the determination of dE/dx of partons in the QGP. 
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Figure 27: Plots from PHENIX [88] of Sp' T /p' T at p' T = pr(p + p) = 7 GeV for 7r° (RHIC) and charged 
hadrons (LHC): a) as a function of centrality (7V part ), b) as a function of dN^/dp. 


8.3 At last: jet measurements in Au+Au at RHIC in 2014? 

Some interesting new jet measurements in A+A collisions at RHIC were presented at Quark 
Matter 2014 in a plenary review talk on jets by Yen-Jie Lee who works on CMS [89]. Figure [28| 
shows that the STAR charged jets in a cone with R = 0.2 have much less suppression 
( Raa ^ 0.3) than 7T° (0.2 < Raa < 0.3) in the range 10 < Pt < 20 GeV. 
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Figure 28: a) (left) STAR Raa for charged jets at NN =200 GeV in central Au+Au collisions (see 
details in legend) compared to b) Raa for PHENIX 7 T °. The dashed line at 0.3 is the maximum Raa for 7 r° 
in this pt range. 


This is quite different from jets at the LHC (Fig. 29) which have comparable or smaller Raa 
than charged particles from jet fragmentation in the range 30 < Pt < 100 GeV. Note that 


the 7, W and Z {) bosons in Fig. [29b which are not coupled to color are not suppressed. 


27 


























Figure 29: a) (left) Raa for jets at ^/s NN =2.76 GeV by CMS and ALICE compared to b) CMS Raa for 
charged hadrons (Raa ~ 0.55), 6-quarks and 3 favorite Electro-Weak Bosons |8F ■ 


For STAR, the disagreement of the jet and single particle Raa gets worse as the jet cone 
is increased from R= 0.2 to 0.3 to 0.4 (Fig. 30). Some people would say that this is great 
because all the jet fragments and/or any energy lost in the QGP by the originating parton 
have been captured in the R=0A cone. Skeptics like myself can hardly wait to see what 
happens when the jet cone is further increased. After 14 runs at RHIC, the jet learning 
curve in Au+Au central collisions still has a way to go. 
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Figure 30: STAR Raa for charged jets at ^/s NN =200 GeV in central Au+Au collisions for 3 different jet 
cones with R = 0.2, 0.3, 0.4 (see details in legend and sketch) [89]. 
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The good news for the future is that a new detector, now called sPHENIX, to find jets by 
the more traditional method using hadron calorimetry has been proposed, is moving along 
on the approval process and is on the schedule at RHIC (Fig. [5]b) for comissioning and full 
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Figure 31: a) (left) new detector in place of PHENIX; b) (right) beam’s eye view with details [90J. 

installation in the PHENIX IR in 2019-20. It is based on the (made in Italy) BABAR 
superconducting solenoid from SLAC which became available when the B-factory in Italy 
was unfortunately cancelled. Design work is moving along quickly (Fig. |3T| so now is the 
optimum time for new collaborators to join. 
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